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1. Introduction 

Entanglement in continuous variables has been of great interest since the celebrated 
paper of Einstein, Podolsky and Rosen (EPR) who constructed a two-particle state 
which was strongly entangled both in position and momentum spaces. The EPR state 
has been the subject of many discussions on the nonlocality of quantum mechanics. 
It turns out that the EPR state can be physically realized in a high gain parametric 
amplifier. This opened up the possibility of a variety of new experiments j21ISllllEl using 
entanglement in continuous variables. The Wigner function for such states is Gaussian 
in position and momentum variables jS]. The Gaussian states are very special in the 
sense that the information on the higher order correlations can be extracted from second 
order correlations. The criteria for entanglement these states has been formulated in 
terms of second order correlations between position and momentum variables [3 
Mancini et al. derived an equivalent set of criteria. A great advantage of these 
inequalities is that the transpose criteria has been translated into something which is 
directly measurable. In this paper we focus our attention on the entangled character of a 
family of non-Gaussian states, viz., the pair coherent states. These states are entangled 
since the expression for pair coherent states is already in Schmidt form. We examine its 
entanglement character in terms of the Peres-Horodecki criteria. We calculate explicitly 
the eigenvalues of the partial transpose of the density matrix and show that some of 
these are negative. We also present results on the correlation entropy and linear entropy. 

To setup an experiment to quantitatively measure entanglement, one can think 
of quasi-probability distribution functions, viz. the Glauber-Sudarshan P-function, 
Q-function, and the Wigner function. We study the relationship between the non- 
classicality the P-function of the state and entanglement. It is known that if P-function 
of the state is well-behaved, then the state is separable. We study the P-function for 
the pair coherent state to show the entangled nature of the state, we also check the 
entanglement in pair coherent state using the inequalities for second order moments. 

The organization of this paper is as follows. In Sec. II, we introduce the family 
of bipartite non-Gaussian states of radiation field. We describe their properties briefly. 
In Sec. Ill, we investigate the inseparability of the pair coherent states in light of 
Peres-Horodecki criterion and von Neumann entropies. Later we discuss the relation 
between entanglement and non-classicality of P-function. Finally we study the existing 
separability inequality to detect entanglement in the pair coherent state. 

2. Pair coherent state: An entangled non-Gaussian state 

The simplest examples of non-Gaussian states of the field are, say, the single photon 
states. Other examples could be states generated by excitations on a Gaussian state 
[Tin ITT] . The state which has been extensively studied for its nonclassical properties 
and violation of Bell inequalities ^211^ is the pair coherent state A pair coherent 
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Figure 1. Contour plot of the quadrature distribution P{xa,Xb) for the pair coherent 
state for = — 1. 



state ICl) is the state of a two- mode radiation field 

ab\Cq) = C\Cq) , 
{a^a-b^b)\C,q)=q\C,q) , 



with the following properties: 

(la) 
{lb) 



where a and b are the annihilation operators associated with two modes, C is a complex 
number, and q is the degeneracy parameter. The pair coherent state for q = 
(corresponding to equal photon number in both the modes) is given by 

oo /-n 



/•r 

\C,0) = NoY.^\n,n) 



n=0 



(2) 



where A^o = l/\/-^o(2|C|) and /o(2|C|) is the modified Bessel function of order zero. The 
coordinate space wave function is given by 



oo 

a,Xb\C,0) = NoY^ —r{xa\n){xb\n) 



n=0 

oo /-n 



No 2^ — — 7^, exp 



n=0 



xl + xl 



(3) 



where {xa\n) is a harmonic oscillator wave function given in terms of the Hermite 
polynomial as 



{xa\n) = (2"n!v^) 



H4xa)e-^^^r^ . (4) 

It is clear from the expression ^ that the wave function of the pair coherent state is 
non-Gaussian. We have shown the quadrature distribution P{xa,Xb) = | (x^, 0) p 
of the state (jHI) in Figure ^ The distribution reflects the entanglement present in the 
state. Note that the pair coherent state can be obtained by projecting the two-mode 
coherent state 



-(|ap + |/3p)/2 



E 

n,m=0 



\n, m) 



(5) 
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onto a space with a fixed difference of the number photons in two modes. The well 
known squeezed vacuum state with Gaussian wave function is given by 



IC)TP = Vl-|CPEC1r^,n). (6) 



ra=0 



Note that the expansion coefficients are different in (j2)) where the coefficients decrease 
quickly with increase in n\. 



3. Inseparability of the Pair Coherent State 

In this section we quantitatively study the entanglement in pair coherent state. Note 
that the state has an obvious form of Schmidt decomposition. This reflects the fact 
that this state is an entangled state. In the next subsections, we examine the other 
criteria to give an estimate of its entanglement. 



3.1. Peres- Horodecki inseparability criteria 

The Peres-Horodecki inseparability criterion |T3] is known to be necessary and sufficient 
for the (2 x 2) and (2 x 3) dimensional states, but to be only sufficient for any higher 
dimensional states. This criterion states that if the partial transpose of a bipartite 
density matrix has at least one negative eigenvalue, then the state becomes inseparable. 
The density matrix p corresponding to the state |C;0) (which is a infinite dimensional 
state) can be written as 

P= (flCnn\n,n)\ C*^^{m,m\\ , (7) 

\n=0 / \m=0 / 

where Cmm. = ^o^- Hence the partial transpose of p is given by 

oo 

PPT = CnnC^Jn,m){m,n\ . (8) 

ra,m=0 

One can now calculate the eigenvalues and eigenf unctions of the matrix pprp as follows: 
Let us start with the following set of two Hermitian conjugate terms for n 7^ m in the 
above equation 

CnnC^m\n,m){m,n\ + CmmC*n\m,n){n,m\ . (9) 
Diagonalizing the above block of the matrix ppT we find the following eigenvalues: 

- /o(2|C|) {n\y ' 
/o(2|C|) n!m! 

and the corresponding eigenf unctions \n,n) and {\n,m) ± e~'^^\m,n)) /\^, where 9 
is the relative phase of the amplitudes C„„ and Cmm and is defined by e*^ = 
CnnCmm/\^nn\\Cmm\- Clearly the matrix pprj has several negative eigenvalues. Hence 
according to the Peres-Horodecki criterion, the pair coherent state is an inseparable 
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Figure 2. Variation of the correlation entropy of the pair coherent state (sohd hne) 
and the squeezed state (dashed hne) with |C|. For large values of \(\ the correlation 
entropy varies linearly with |^|. 



state§. Note that if the phase of the parameter ( is random, then the state becomes 
separable, as then terms corresponding to different values of n and m drop out of the 
double summation in (jH)). 



3.2. Correlation entropy 

Correlation entropy of a bipartite system consisting of subsystems a and b is given by 

m 

-^corr = Sa -\- Sh — Sab ; (H) 

where S^ is the von Neumann entropy of the system k. If a and b are uncorrelated 
(separable), then /corr vanishes. Now for any bipartite pure stcitc, Sfii) is zero. We have 
calculated Sa^b for the pair coherent state as 

c^-c,- ^ ICP" , ^ f KP" ^ .... 

- - h mc\)^' wmW") ■ ^ ^ 

We plot the correlation entropy ()11|) for the pair coherent state with \(\ in Figure |21 
The correlation entropy for the pair coherent state remains non-zero for all values of |C| 
which signifies that the state is inseparable (entangled) for all \(\. For smaller values 
of Id, the entropy increases slowly; but at larger values of it saturates. Note that 
the squeezed vacuum state (jHl) is a very special kind of Gaussian non-classical state, 
entanglement properties of which have been much studied in literature j^]. For a better 
understanding of the inseparability of the pair coherent state, we have compared the 
correlation entropies of the pair coherent state and of the squeezed vacuum state in 
Figure El Clearly the correlation entropy of the squeezed vacuum state grows much 

§ The eigenvalues of the partial transpose of the density matrix for the squeezed vacuum state are also 
given by with nl and ml replaced by unity. 
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Figure 3. Variation of the linear entropy Iim for pair coherent state (solid line) and 
squeezed vacuum state (dashed line) with 



faster than that of the pair coherent state. This is because, as |C| approaches to unity, 
the squeezed vacuum state becomes much more incoherent than the pair coherent state. 
At Id = 1 the squeezed vacuum state becomes completely random, as all the possible 
states of the either mode then become equally probable. It is well-known that the 
correlation entropy of a system with all basis states equally probable is 2 log2 A^, where 
is the number of possible basis states of the system. For squeezed vacuum state, 
as A^ ^ oo, the correlation entropy diverges for \(\ = 1. On the other hand, in case 
of pair coherent state, the picture is different at |C| = 1. In this case the states with 
lower occupation number become more probable than the states with higher occupation 
number. Thus the correlation entropy for the pair coherent state remains less than that 
of the squeezed vacuum state for 1^ approaching unity, suggesting that the pair coherent 
state remains much more coherent than the squeezed vacuum. 



3.3. Linear entropy 

We further calculate the linear entropy of the pair coherent state, which is given by 
/lin = 1 — Tr(p|), where pk is the reduced density matrix of the subsystem k. For a 
pure state density matrix p, In^ vanishes as Tr(p^) = 1. But for an entangled state, pk 
does not have the form of a pure state density matrix. Thus, any non-zero /un provides 
signature of entanglement present in the state. Note further that the linear entropy is 
closely related to the entanglement measure in terms of Schmidt number For a 

pair coherent state, the linear entropy is given by 

oo I ^ 1 4n 



'lin 



/o(2|C|)^„" 



E ^ • (13) 



=0 



We show the variation of the quantity Jun with \(\ in Figure IHl Clearly for any non- 
zero Id, /lin becomes non-zero and positive implying entanglement in the pair coherent 
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Figure 4. Variation of product of the variances of the joint position and momentum 
of two subsystems in a pair coherent state (sohd hne) and in a two-mode squeezed 
vacuum state (dashed Hne) with the squeezing parameter |C|, for (p = n and m — I. 
The product remains less than unity for all non-zero |C|- 
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Figure 5. Variation of total variance in the EPR-like variables u and v with |C| for 
pair coherent state (solid line) and squeezed vacuum state (dashed line) for m = 1. 
We have chosen (j> = ir. 

state. For |^| ~ 0.4, Jun ~ 0.25, i.e., Tr(p|) is close to unity. This implies that the state 
represented by pk is more like a pure state than a mixed state. On the other hand, 
for ICI = 1, Jiin ~ 0.6, which means that pk represents more like a mixed state than a 
pure state. Thus the state Q is more entangled. Further we have plotted the linear 
entropy of the squeezed vacuum state © in Figure El Clearly for small \(\ the linear 
entropy of both the states show similar behavior. But at larger values of \(\, linear 
entropy of squeezed vacuum state is more than the pair coherent state, i.e., the degree 
of mixedness of the pair coherent state remains lower than the squeezed vacuum state 
for larger values of \(\. 
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3.4- Entanglement and non-classicality of the P function 

So far we have investigated the entanglement in pair coherent state in terms of Peres- 
Horodecki inseparabihty criterion and various entropies. However all these criteria are 
not possible to verify in experiments. We will now the study the non-classicality of P- 
function for the pair-coherent state as P-function of a density matrix can be measured 
in experiments. For bipartite systems, the two-mode density matrix can be written in 
terms of the diagonal coherent state representation as 

p = J J d^ad'^(3P{a,a*;f3,(3*)\a,(3){a,(3\ . (14) 

It is known that if the P-function has non-classical character, then the state is entangled. 
We examine the inseparability of the pair-coherent state from the point of view of the 
non-classicality of the P-function. The Glauber- Sudarshan P-distribution function gives 
a quasi-probability distribution in phase space, which can assume negative and singular 
values for non-classical fields. There is another distribution function called Q-function 
which is related to the P-function by 

Q{a,a*] P, P*) = -^{a, /3\p\a, (3) 

= ^J P(7, 7*; 6, 5*)e-\''-^\\~\'-^\"d^^dH , (15) 

which is always positive. Note that if the function P(7; 5) were like classical probability 
distribution, then Q(a; /3) > V However if Q is zero, then P must become at 

least negative (referring to nonclassicality of the state) in some parts. Hence the exact 
zeroes of the Q-function are also a signature for the non-classicality of the field. In order 
to see the non-classicality of the pair coherent state, we examine the structure of the Q 
function. The Q function for the pair coherent state can be calculated as 

Q(a, «*;/?,/?*) = ^^e-(l"l'+l^l')|Jo(2yC^)|^ . (16) 

This function has zeroes only if a and (3 are out of phase for real positive and if 
2^^|a||/5| = zq, where zq are the exact zeroes of the Bessel function Jq{z). The smallest 
few values of zq are 2.4048, 5.52, 8.6537, 11.7915, 14.9309 etc. Existence of these zeroes 
proves that the pair-coherent state is a non-classical state. 

However it is worth mentioning that, the notion that, if the P-function of a state 
is not well-behaved, then the state in inseparable, is not always true. For example, for 
a two- mode separable Fock state \n,m), the P-function is not at all well-behaved. This 
motivates us to investigate some alternative inseparability criteria, which can be verified 
in experiments. We will show that existing inseparability inequalities are quite useful 
in detecting entanglement in non-Gaussian state like in experiments, albeit under 
certain condition. 

3.5. Separability inequalities 

Duan et al. |7^ and Simon jSj independently have derived the separability criterion of 
a bipartite continuous-variable system in terms of the second-order correlations. This 
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criterion states that if a state is separable, then the uncertainties in a pair of EPR-hke 
operators u and v satisfy, 

M = {(Auf) + {(Avf) >m' + ^, (17) 

where 

u = \m\Xa H Xf, , (18a) 

m 

V = \m\pa Pb , (18&) 

m 

for any arbitrary nonzero real number m. Here Xk = {k + k"^) / \pl and Pk = (k — k^)/iy/2 
{k = a, b) are the position and momentum operators for the subsystem k satisfying the 
commutation relation [xk,Pk'] = i^kk'- Conversely, violation of this criterion provides a 
sufficient condition for inseparability of states, albeit with a lower bound 

< M < + (19) 



^ 



9 ' 



which for m = ±1, reads as 

< M < 2 . (20) 

For a bipartite Gaussian state, the criterion (fTTjl is also sufficient for separability. 

Equivalent necessary and sufficient condition for separability of Gaussian states 
have been derived by Englert and Wodkiewicz ^H] using density operator formalism. 
They have shown that the positivity of the partial transposition and P-represent ability 
of the separable Gaussian states are closely related. These criteria have been 
experimentally verified via the interaction of linearly polarized field with cold atoms 
[in], in atomic ensembles [211!, and with squeezed hght fields |2I1I221- Mancini et al. [Hj 
have shown that separability of a state leads to the following uncertainties in a pair of 
EPR-like variables 

M. = ((An)2)((At;)2)>l, m = l, (21) 

where u = Xa + xi, and v = Pa — Pb- Violation of this inequality provides a sufficient 
criterion of inseparability in Gaussian states. 

We will now discuss the validity of the criterion ()21|) in case of a pair coherent state 
which is a non-Gaussian state. We calculate the uncertainties ((Au)^) and ((Af)^) (for 
m = 1) for the pair coherent state |C, 0). We find the averages {xa + x^) = and 

{{xa + Xbf) = (C|(l + a^a + b^b + ab + a^b^)\0 

= l + 2|Cim] + 2|C|cos0. (22) 

Here Ii{2Q is the modified Bessel function of order one and ( = ICIe**^. Thus 
{{AuY) = {v?) — (u)^ = {{xa + XbY). In a similar way, one can calculate the variance 
((Af )^) which is found to be equal to ((Am)^). In Figure|31 we have shown the variation 
of product {{Auf) {{AvY) with Clearly, the inequality ((211 is violated for all 
We can thus infer that the pair coherent state is an inseparable state. 
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Now we will discuss whether the criterion [Equation (jl9|) ] is applicable for the pair 
coherent state. The total variance M = ((Am)^) + ((Af)^) can be calculated as 



M= \m\' + — )+2[ |m|' + — Id 



mV V ' mV' 7o(2|C|) 



+ 4^|C|cos0. (23) 



m 



We show the variation of the above quantity with |C| in Figure El for = tt. The figure 
shows that total variance remains less than |mp + ^ for all \(\. Thus the inequality 
(fTTHl is satisfied for the pair coherent state under the condition 

(sign of m)(sign of coscp) < . (24) 

Thus the criterion ()19|) is sufficient for inseparability of the pair-coherent state only if 
the above condition is satisfied. 

We have also compared the degree of violation of (PT|) and ()19|) of the pair-coherent 
state with that of the squeezed vacuum state (jH)) in Figure IH and Figure El It shows 
that the criteria ^IT\ and ()19|1 are violated more by the squeezed vacuum state for larger 
\(\. Thus the degree of inseparability is more for the squeezed vacuum state at larger 
\(\. This is due to the fact that the expansion coefficients in the pair coherent state [see 
Equation decrease quickly with the increase of n\. 



4. Conclusions 



In conclusions, we have studied the inseparability of a special family of non-classical 
states, called the pair coherent states which are non-Gaussian in nature. We confirmed 
the inseparability of pair coherent states in the light of Peres-Horodecki criteria and 
various entropies. We then demonstrated that the existing inseparability criterion ()19j) 
based on second order correlation, is applicable to these kind of non-Gaussian states 
only under certain constraints. 
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